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SunTO. Si considera il problema di Cauchy per sìstemi strettamente iperbolici di leggi di conservazione 
us + [F(u)], =0, u(0,7) = U(2), (PC) 


in una singola variabile spaziale r € IR Si discute la teoria generale sull’esistenza di soluzioni di (PC), 
globalmente definite nel tempo, per dati iniziali 7 = U(x) con variazione totale sufficientemente piccola. In 
particolare, si presentano alcuni risultati recentemente ottenuti in collaborazione col Dott. Andrea Marson, 
per una classe di sistemi che non soddisfano classiche ipotesi di convessità della funzione F = F(u). 


ABSTRACT. We consider the Cauchy problem for strictly hyperbolic systems of conservation laws in one 
space dimension 
u + [F(u)], =0, u(0,7) = U(z). (CP) 


We discuss the general theory on the existence of global solutions to (CP), with initial data with small total 
variation. In particular, we present some recent results obtained in collaboration with Dott. Andrea Marson 
for a class of systems which do not satisfy classical convexity assumptions on the flux funxtion F = F (u). 
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1. INTRODUZIONE 


In questo seminario consideriamo il problema di Cauchy per sistemi iperbolici di leggi di 
conservazione 


u + [F(u)], =0, (1) 
u(0,2) = (2), (2) 


in una singola variabile spaziale r € R. Qui u= u(t,r) può essere interpretata come la densità di 
una quantità fisica conservata e F = F(u) come la corrispondente funzione flusso. 

Discuteremo la teoria generale sull'esistenza di soluzioni di (1)-(2), globalmente definite nel 
tempo, per dati iniziali © = (1) con variazione totale sufficientemente piccola. In particolare 
vogliamo presentare alcuni risultati, recentemente ottenuti in collaborazione col Dott. Andrea 
Marson, per una classe di sistemi che non soddisfano classiche ipotesi di convessità della funzione 
F = F(u), nel senso che verrà meglio precisato in seguito. 


La funzione flusso F : 2 + RR" si suppone essere regolare e definita in un intorno dell’origine 
9 Cc PR". Denotiamo Ai(u),...,An(u), gli autovalori della matrice Jacobiana A(u) = DF(u), e 
chiamiamo ri(u),...,fn(u), una base di autovettori destri di Alu). 

L'ipotesi generale che faremo su (1) è che il sistema sia strettamente iperbolico e cioè che 
gli autovalori A;(u) siano tutti reali e distinti. In particolare possiamo supporre che 


Ma (U) < --- < An(u) Vue. (3) 


Una caratteristica fondamentale di equazioni nonlineari del tipo (1) è che, in generale, la 
soluzione del problema di Cauchy sviluppa delle discontinuità in tempo finito, anche nei casi in cui 
si considerino condizioni iniziali regolari (di classe C°°). Pertanto, volendo ottenere dei risultati di 
esistenza globale ed essendo il sistema (1) scritto in forma di divergenza, risulta naturale lavorare 
nella classe delle funzioni discontinue interpretando la soluzione di (1) come soluzione nel senso 
delle distribuzioni. Più precisamente, diremo che una funzione localmente integrabile u = u(t, x) è 
soluzione debole del problema di Cauchy (1)-(2) se 


n (n fu - © + F(u)- $2]drdt + (C T(2) -$(0,2) dr =0 4 


per ogni $ € C* a supporto compatto. Consideriamo ora una soluzione debole u = u(t, 7) di (1.1) 
che sia discontinua lungo una retta (nel piano z-t) di equazione £ = r(t) e chiamiamo 


S(t)= li t, D(i)= li t,2), 5 
u°(t) a 2), u (8) a 2) (5) 
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i limiti sinistro e destro di u(t, -) nel punto di salto. Fissando un’arbitraria funzione test $ € C! con 
supporto in una regione I' attorno a 7, e applicando il teorema della divergenza al campo vettoriale 
(u- $, F(u) - $), si può facilmente dedurre dalla (4) che lungo 7 deve valere la relazione 


(uP — u5) 4= F(uP)- F(u5). (6) 


L'equazione vettoriale (6) è nota come condizione di Rankine-Hugoniot (RH). Data una 
coppia di stati (u5, uP), diremo che uS, uP soddisfano le condizioni (RH) se verificano la (6) 
per un qualche scalare A = 4 € RR. Fissato uno stato (sinistro) u° € O, per determinare l'insieme 
degli stati (destri) u = (u1,...,tn) € £ che soddisfano la (6) dobbiamo risolvere il sistema di n 
equazioni scalari 


(u-u?)\= F(u)- F(u°) (7) 


nelle n + 1 incognite (u1,...,un,A). Usando il teorema delle funzioni implicite, si mostra che 
(localmente) tale insieme consiste di n curve regolari. Più precisamente, esistono n curve regolari 
o + S;(u®)[0], i=1,...,n, ed n funzioni scalari o + \î(u®)[o], i= 1,...,n, che soddisfano 


(S;(0°)[0] - u°) A (u9)[o] = F(5;(u°)[0]) — F(u°) Vo. (8) 


La curva o + S;(u®)[0] si chiama curva di shock dell’i-sima famiglia (passante per u°), mentre 
Aî (®)[o] sarà chiamata velocità dello shock (u°, S;(u°)[c]). Si noti che il vettore S;(u°)[o]-u® 
risulta essere un’autovettore della matrice media 


1 
Au, u”)) = ti A(u” + é(u” — u°))de, u” = S;(u®)[g], 


avente come autovalore corrispondente la velocità dello shock A#(u°)[c]. Se ne deduce che la velocità 
dello shock assume un valore vicino alla media tra le velocità caratteristiche Ai(u®) e A;(u”). In 
particolare, si ha il seguente sviluppo 


A (u9)[0] = A;(u®) + 3DA:(u9) -ri(u°) + O(1) [of?, 


dove il simbolo di landau ©(1) denota una quantità uniformemente limitata in valore assoluto da 
una costante indipendente dal sistema (1). 


Possiamo ricapitolare le considerazioni fin qui fatte dicendo che, se u = u(t,r) è una soluzione 
debole di (1), allora in ogni suo punto di discontinuità il corrispondente stato destro uP (definito 
come in (5)) deve appartenere ad una delle n curve di shock passanti per lo stato sinistro uS, e la 
velocità della discontinuità deve essere precisamente la velocità dello shock. 


In realtà le condizioni di Rankine-Hugoniot sono in generale insufficienti per determinare 
un'unica soluzione debole del Problema di Cauchy (1)-(2), come mostra il ‘seguente 
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Esempio 1. Il problema di Cauchy per l'equazione (scalare) di Burger 


uz 
Ut + (E. = 0, (9) 
0 se t<0, 
0,r)= 10 
UO) { 1 se z>0, La) 


ammette infinite soluzioni deboli. Precisamente, per ogni fissato a € [0, 1), la funzione 


0 se r<at/2, 
u°(t,r) = a se  at/2<x<(1+0)t/2, (11) 
1 se r>(1+0)t/2, 


x=0.t/2 


x= (0+1)t/2 


0 


FIGURA 1 


fornisce una soluzione di (9)-(10) nel senso delle distribuzioni. Infatti ua è costante al di fuori delle 
due linee di discontinuità €&2(t) = at/2, €9(t) = (1+a) t/2, e dunque soddisfa l’equazione (9) nelle 
regioni non contenenti le linee £f, 7. D'altra parte si verifica immediatamente che ua soddisfa 
anche le condizioni di (RH) lungo le due discontinuità £7, €. 


Pertanto, allo scopo di determinare un’unica soluzione del problema di Cauchy (1)-(2), è 
necessario aggiungere al concetto di soluzione debole delle ulteriori condizioni di ammissibilità, 
possibilmente motivate da considerazioni fisiche, che permettano di selezionare come ammissibili 
solo certi tipi di discontinuità. Varie condizioni di ammissibilità sono state proposte in letteratura 
dando luogo a differenti definizioni di soluzione: 

1) Limite viscoso: la soluzione è esprimibile come limite di soluzioni di equazioni paraboliche 

(viscose) del tipo 

uf + A(u°)uf = E USz) 

al tendere a zero del termine viscoso; 

2) Disuguaglianza entropica: la soluzione soddisfa nel senso delle distribuzioni una disuguaglianza 

associata ad una coppia di funzioni (7, g), in cui 7: R° + R è un’entropia convessa per (1) 

e q:R" + R èla corrispondente funzione flusso entropico; 
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3) Condizione di stabilità di Laz: in ogni punto di discontinuità (uS, uP) della soluzione, se gli 
stati u e uP sono connessi da uno shock della i-sima famiglia che si muove con velocità s, 
allora lo shock è ammissibile se si ha 


Ai(u5) > 5> Ai(u!). (12) 


Nel prosieguo della discussione, ci limiteremo a considerare solo quest’ultimo tipo di condizioni di 
ammisibilità e dunque chiameremo soluzione entropicamente ammissibile (alla Lax) di (1) 
una soluzione debole di (1) che, in ogni punto di discontinuità, soddisfi le condizioni di Lax. Si noti 
che, nel caso dell’esempio sopra considerato, nessuno dei due shock €f, €, soddisfa le condizioni 
di Lax (12), per qualunque a € [0, 1]. 


FIGURA 2 


2. IL PROBLEMA DI RIEMANN 


Tre metodi si sono sviluppati in letteratura per costruire soluzioni deboli entropiche di (1)-(2) 
globalmente definite nel tempo: 
i) Lo schema di Glimm [G], [L13]. 
ii) Il tracciamento dei fronti d'onda [B1], [DP1], [Ri]. 
ili) La viscosità evanescente [DP2]. 


Questi metodi si basano tutti su argomenti di compattezza: stime a-priori della variazione to- 
tale e un’applicazione del teorema di Helly (in i) e ii)) o risultati di compattezza compensata 
(in iii)) garantiscono che qualche sottosuccessione di soluzioni approssimate converge in L},, ad 
una soluzione debole, globalmente definita per ogni tempo t > 0. Noi concentreremo la nostra 
attenzione sul metodo del “tracciamento dei fronti d’onda” (wave front-tracing) che, assieme allo 
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schema di Glimm, fornisce i risultati più generali essendo il metodo della “viscosità evanescente” 
applicabile solo per.sistemi di n < 2 leggi di conservazione. 


Il mattone principale per la costruzione di soluzioni approssimate, sia nel tracciamento dei 
fronti d’onda che nello schema di Glimm, è costituito dalla soluzione del problema di Riemann, 
cioè.del problema di Cauchy con dato iniziale © avente un'unica discontinuità nell’origine, e valore 
costante a sinista e a destra dell'origine: 


uS se T<0 SD 
u(c) = ’ O. 1 
ha { uP se r>0, iaia sa 


La soluzione u = u(t, x) di questo tipo particolare di problema di Cauchy si esprime come funzione 
auto-simile, i.e. u(t,z) = w(7/t) per una qualche.w: R+ R", eventualmente discontinua. Prima 
di discutere l’espressione generale della soluzione di (1)-(13), dobbiamo introdurre delle ulteriori 
ipotesi sugli autovalori \;(u) e sugli autovettori ri (u) della matrice Jacobiana DF(u) della funzione 


flusso. A tale scopo diamo la seguente 


Definizione 1. L’i-simo campo caratteristico rj(u) si dice 


(GNL) genuinamente non-lineare se 
DA;(u)-ri(u) #0 VueD; 
(LD) linearmente degenere se 
DA;(u)-ri(u)=0 VueD 


(NGNL) non genuinamente non-lineare (con un solo punto di fiesso) se non è (GNL) 
e l’insieme 


r= {ue 0 : DA(u)-ri(u)=0} (14) 


è una varietà n — 1 dimensionale di Q, trasversale al campo caratteristico r;, tale che 


D(DA;-ri)(u)-ri(u) #0 Vuerì. (15) 


Si noti che le ipotesi (LD) o (GNL), nel caso in cui n= 1 (cioè per equazioni scalari), cor- 
rispondono a richiedere, rispettivamente, che la funzione flusso F = F(u) sia affine o strettamente 
convessa (o concava). Mentre l’ipotesi (NGNL) equivale a richiedere che F = F(u) non sia nè 
concava, nè convessa, ma abbia un-unico punto di flesso. Per semplicità, ci riferiremo dunque a 
quest’ultima ipotesi come all'ipotesi di non convessità della funzione flusso. ‘ 


Per il prosieguo della discussione, è utile introdurre anche la seguente notazione relativa ad i 
campi caratteristici ri,î = 1,...,n. Fissato u® € Q, chiameremo curva di rarefazione dell’i- 
sima famiglia (passante per u°) la curva integrale di ri attraverso u°, i.e. la curva ottenuta 
risolvendo il problema di Cauchy nello spazio degli stati 


du 
rta ri(u), u(0)=u0. 
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Denoteremo tale curva come o + R;(u)[0]. Una proprietà fondamentale che lega le curve di shock 
e di rarefazione della stessa famiglia, passanti per lo stesso punto u°, è che esse risultano avere 


(1) 


in u° un contatto del second’ordine. Più precisamente, usando un’opportuna parametrizzazione, 


si ottengono i seguenti sviluppi 


Ri(u°)[o] = +ori(u0) + °° Dr;(u9) -ri(u°) + O(1)[o}}, 
3 (16) 
S;(u°)[0] =u + ori(u°) + FDri(u®) . ri(u0) + O(1)[o]}. 


Consideriamo ora i diversi tipi di onde elementari che compongono la soluzione generale 
di (1)-(13) a seconda delle diverse ipotesi fatte sui campi caratteristici. 


1. Onde centrate di rarefazione. Supponiamo che l’i-simo campo caratteristico sia (GNL) 
e, in particolare, possiamo scegliere l'orientamento di r; in modo tale che valga 


DA;(u)-ri(u)> 0 VueEA. (17) 


LIA uP € Q appartenente alla parte positiva della curva integrale di r; passante per 
SEG, ie. 
PD = Ri(u5)[o], per qualche o >0. (18) 


Per ogni a € [0, 0], definiamo 
Ai(2) = A;(R;(u°)[q]). 
A causa di (17) la funzione a + \;(a) è strettamente crescente. Possiamo allora definire la funzione 


continua, regolare a tratti 


uS se r<t);(uS), 
u(t,r)= Ri(u5)[a] se r=t);(0), a € [0, o], (19) 
uP se T>tA;(uP), 


che risulta essere una soluzione debole entropica del problema di Riemann (1)-(13). Infatti 
l'equazione (1) è ovviamente soddisfatta nelle regioni in cui z < tA;(uS)ox>tA;(uP), dove u è 
costante e dunque w; = uz = 0. D'altra parte, poichè u è costante lungo ogni retta di equazione 
c=tA;(a), @ € [0, o], si ha 


us(t,2) +A;(a)uz(t,7) =0 z=t);(a). (20) 


Osservando allora che 


du(t,2) _ dR;(uS)a] (dAi(a)\-1_dA;(t,2) dAi(a)\-1 1 
da" di (da) e e) 


è un autovettore della matrice A(u(t,z)), con autovalore A;(a) = A; boa dalla (20) segue che 
u è soluzione dell'equazione (1) scritta in forma quasi-lineare 


u + A(u)uz = 0. (21) 
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Si noti che l’ipotesi o > 0 in (18) è fondamentale per garantire la buona definizione di u in quanto, 
nel caso o < 0, la funzione definita da (19) risulterebbe essere multivoca nella regione in cui 
c/t e [A;(uP), A;(uS)] (assume in ogni punto tre valori). 


2. Shock. Supponiamo ancora che l’i-simo campo caratteristico sia (GNL) e che valga 


la (17). Consideriamo uP € Q appartenente alla parte negativa della curva di shock della i-sima 
famiglia passante per u5 € Q, i.e. 


uP = S;(u5)[0], per qualche o <0. (22) 


Allora, chiamando A° = Af(u5)[0] la velocità dello shock determinata dalle condizioni di (RH) (8), 


definiamo la funzione é 


vu se T<KEA 
u(t,2)= { DI i 


23 
u se DISTA, (28) 


Si noti che, essendo o < 0 e poichè vale (17), la velocità caratteristica A; è monotona decrescente 
lungo la parte negativa della curva di shock e dunque si ha 


Ai(uP) < A8(uS)[a] < Ai(u5). (24) 


Ne consegue che la funzione definita da (23) è una soluzione debole entropica nel senso di Lax. Nel 


caso in cui uP = S;(u5)[c] con o > 0 invece, la funzione definita da (23) non soddisfa le condizioni 
di stabilità di Lax. 


3. Discontinuità di contatto bilatere. Supponiamo che l’i-simo campo caratteristico 
sia (LD). Consideriamo uP € Q appartenente alla curva integrale di r; passante per uS € D, i.e. 


uP = R;(u5)[0], per qualche o. (25) 


Per ipotesi la i-sima velocità caratteristica A; è costante lungo questa curva. Ponendo allora 
A° = A;(uS), la funzione definita da (23) risulta essere una soluzione debole del problema di 
Riemann. Infatti, in questo caso si ha 


F(uP)- F(u5)= : DF(R;(u5)[a]) -r;(Ri(uS)[a]) da 
"30 (26) 
=). (Ri(u9)o] - uS), 


da cui si deduce che u5, uP soddisfano le (RH), qualunque sia il valore di o in (25). D'altra parte 
le condizioni di Lax sono banalmente soddisfatte in quanto si ha 


Ai(uS) = 1° = Ai(uP). (27) 


Si noti che la relazione (26), in particolare, mostra che le curve di shock e di rarefazione coincidono 
per campi caratteristici (LD). 
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4. Onde miste: discontinuità di contatto unilatere + onde centrate di rarefazione. 
Supponiamo che l’i-simo campo caratteristico sia (NGNL) con un solo punto di flesso e, in par- 
ticolare, supponiamo che valga 


D(DA; * ri) (u) . ri(u) <0 Vuer,, (28) 


cioè che, per ogni fissato w € T;, la funzione a + X; (Ri(u)[a]) assuma il suo massimo in a = 0. 
Consideriamo due stati uS, uP € Q situati da parte opposta rispetto alla varietà F'; definita in (14), 
e supponiamo che tra uS e uP ci sia un punto u* con le seguenti proprietà: 

i) uM appartiene alla parte positiva della i-sima curva di rarefazione per ul € Q, i.e. 


uM = Ri(uP)[01] per qualche 01 > 0; (29) 
ii) uP appartiene alla parte positiva della i-sima curva di shock passante per uM € Q, i.e. 

uP = S;(uM )[c2] per qualche 02 > 0; (30) 
iii) la velocità dello shock (uM, uP) è uguale alla velocità caratteristica Vi(uM), i.e. 


Ai(8) = Ag (uM)[02]. (31) 


Dalle considerazioni fatte ai punti 1. e 2. segue che, ponendo 
Ai(a) = A;(R;(u5)[a]) a € [0, 01), 


la funzione 


uS se r<t);(uS), 
u(t,x)= | Ri(u5)[a] se r=tA;(0),  @E€[0, 01], (32) 
uP se c>t\;(uM), 


risulta essere una soluzione debole entropica costituita da un'onda di rarefazione centrata, attaccata 
ad una discontinuità di contatto unilatera (sinistra). La discontinuità di contatto (uM, uP) soddisfa 
le condizioni di Lax 

di(uM) = A7(u*)[02] > Ai(u?). (33) 


Fissato uno stato u® € I, per determinare l'insieme degli stati (sinistri e destri) uS, uP € 0, uS € 
R;(u®), che si trovano da parti opposte rispetto a I, e per cui esista uno stato intermedio ul con 
proprietà del tipo (29)-(31), dobbiamo risolvere il sistema di n equazioni scalari 


(u- R:(9)[a])A:(Ri(u9)[0]) = F(u) -— F(R;(u°)(a]) (34) 


nell’incognita (a, u). La relazione (34) implica che il punto u si trova sulla curva di shock che 
passa attraverso il punto R;(u°)[a], e che la velocità dello shock che connette R;(u°)[a] con u è 
precisamente A; (R;(u°)[a]). Pertanto l'equazione vettoriale (34) è equivalente all’equazione scalare 


A3(R:(u°)fa]) [i] = A:(R:(°)la]) (35) 
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nell’incognita (a, t). Infatti (a, t) è una soluzione di (35) se e solo se (a, Si(R;(u°)[a])[t}) è 
soluzione di (34). Chiaramente una soluzione banale di (35) è (a, t) = (0, 0). Usando il Teorema 
della perturbazione singolare di Newton-Kantorovich [D], si dimostra l’esistenza di una soluzione 
non nulla di (35). Più precisamente, si dimostra l’esistenza di una funzione scalare o + pi(u°)[0], 
tale che 


x (R:(00)fo — po] ) [so] = Ai) fe — s(")1]). ‘ (36) 


La curva regolare o + M;(u®)[o] definita da 


Mi(u9)[o] = Si (Ri(u9)fo — us(u)0]) [i (9){0], (37) 


si chiama curva mista dell’i-sima famiglia (passante per u®). Questa curva, come la curva 
di shock o + S;(u9)[0], ha un contatto del second’ordine in vo con la curva di rarefazione 
o + R;(u°)[o]. Una costruzione del tutto simile viene fatta nel caso in cui il campo caratter- 
istico r; sia (NGNL) (con un solo punto di flesso) e valga la 


D(DA; . ri) (u) . ri(u) > 0 VuerT,, (38) 
al posto della (28). 


Introduciamo adesso, per ogni famiglia caratteristica, delle curve definite in termini delle curve 
di shock, delle curve dirarefazione e delle curve miste, che saranno usate per costruire la soluzione 
generale del problema di Riemann. Per ogni fissato uf € O, i € {1,...,n}, 
a) se il campo caratteristico r; è (GNL) o (LD), definiamo la curva 


Ri(u5)o) se  0>O, 
S;(uS)[0] se o <0; 


v:(9)0] = { (39) 


b) se il campo caratteristico r; è (NGNL) (con un solo punto di flesso), soddisfa la (28), e u5 si 
trova a sinistra di T;, allora definiamo la curva 


S;(uS)[0] se o>v(u5), 
PETE Mi(mi(u5))[o+d;(u5)] se — d;(u5) < a < vi(u5), 
tizio R;(u5)[0] se 0<0<-g(u"), dai 
S;(u5)[0] se 0o<0, 


dove r;(uS) = F; N R;(uS), e d;(u5) rappresenta la distanza con segno, lungo la curva di 


5 e mi(uS), i.e. 


rarefazione, tra u 
Ri(u5)[6;(u5)] = ri(uî) e C;. (41) 


Una definizione del tutto simile di V;(u5) viene data nel caso in cui uS si trovi a destra di I, 
o nel caso in cui valga la (38) al posto della (28). 
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La curva o + W;(uS)[o] definita da (38) o da (39) si chiama curva elementare di Riemann 
dell’i-sima famiglia (passante per uS). 


Pertanto, dati due stati uS, uP e 9, se uP appartiene all’i-sima curva elementare di Rie- 
mann Y;(u5), la soluzione del problema di Riemann (1)-(13) risulta essere un'onda elementare: 
una rarefazione, uno shock, una discontinuità di contatto (bilatera), o un’onda mista, a seconda del 
tipo di campo caratteristico. Se invece uP non appartiene ad alcuna curva elementare di Riemann, 
possiamo trovare degli stati intermedi wy = uS, Wi, «.., n = P, in modo tale che ogni coppia 


di stati adiacenti wi_1, w; sia connessa da un'onda elementare, i.e. 


Wi Vi(wi-1)[0;], 1= 1, BSLOA (42) 
Questo può essere fatto per u5, uP , sufficientemente vicini, poichè in tal caso un’applicazione del 
teorema delle funzioni implicite garantisce l’esistenza di onde T1,-:-:0n, univocamente determi- 
nate, tali che 5 
uP = Un (00 ( siro o) [on] . (43) 
D= 
3 u 03 
(07) 
02 
si 
e ga a 
u=%0 1=P1(00)[01] 
FIGURA 3 


A sua volta la (43) determina gli stati intermedi w; che soddisfano la (42). La soluzione generale 
di (1)-(13) può essere dunque ottenuta sovrapponendo le soluzioni degli n problemi di Riemann 


ur + [F(u)], =0, 


0 Wi-1 se x<0, (44) 
ica Wi se c>0, 


su settori differenti del piano x-t. Per costruzione, ogni singolo problema di Riemann ammette una 
soluzione debole entropica costituita da un’onda elementare dell’i-sima famiglia che viaggia con la 
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velocità caratteristica 
1. 


Ae € DET, Ap], AT =Ai(wi-1), Xft= (wi), 


se la soluzione di (1)-(44) consiste di un'onda centrata di rarefazione; 


x = AT Magra 


se la soluzione di (1)-(44) consiste di uno shock; 


A € DET, A], Ne =Awia) ATEA(), 


se la soluzione di (1)-(44) consiste di un’onda mista costituita di un’onda centrata di rare- 
fazione che connette w;-1 ad uno stato W%; (situato tra w;-1 e W;), seguita da una discontinuità 
di contatto unilatera (sinistra). 


Si noti che, supponendo 01,...;0n; sufficientemente piccoli, se ne deduce che le velocità caratter- 
istiche AY delle singole onde elementari risultano essere vicine agli autovalori Aj(wi-1). Pertanto, 
la stretta iperbolicità (3) e la regolarità di F(u) garantiscono che gli intervalli De, A#*4] sono 
disgiunti, i.e., 


ve wWt n= st ,) st 
ae a AT OM x AL, 


. (è ni 
x= La 


FIGURA 4 


Allora, la soluzione (debole entropica) del problema di Riemann ottenuta sovrapponendo le soluzioni 
dei singoli problemi di Riemann (1)-(44), è ben definita. Ad esempio, se i campi caratteristici 
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sono (GNL) e, in particolare, l’i-sima onda è una rarefazione, essa assume l’espressione 


S 


u5 = wo se E<IATT, 
tane Ri(wi-1)(o] se r=tA;(0), a € [0, ci], (45) 
? Wi se ix Sv <A: 
uP=wn se T> LAM. 
Si noti che questa soluzione è una funzione auto-simile, avendo la forma u(t,2) = v(c/t), 


con w : RR" + RR" in genere discontinua. 


3. IL PROBLEMA DI CAUCHY 


Descriveremo adesso l’idea generale del metodo del tracciamento dei fronti d’onda (wave 
front-tracking) usato per costruire soluzioni approssimate di (1)-(2). Una soluzione approssimata 
ottenuta con tale metodo è una funzione u = u(t, x) costante a tratti, che ammette discontinuità 
lungo un numero finito di segmenti o semirette (che chiameremo fronti d’onda) nel piano z-t, e 
che soddisfa nei punti di discontinuità delle condizioni di (RH) approssimate. 


FIGURA 5 


La costruzione inizia al tempo t = 0 prendendo una funzione costante a tratti u(0, -) che ap- 
prossima (nella norma L' ) il dato iniziale a variazione totale piccola. Siano r1 < --- < cy ipunti 
in cui u(0, -) è discontinua. Per ogni £ = 1,..., N, e per tempi positivi in un intorno di (0, ze), ri- 
solviamo in maniera approssimata il problema di Riemann generato dal salto (u(0,z2—), u(0, Ze+)) 
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con una funzione della forma u(t,7) = v((x — xe)/t), per una qualche 4 : R° + R” costante 
a tratti. Più precisamente, se la soluzione esatta del problema di Riemann contiene solo shock o 
discontinuità di contatto, allora u coinciderà con la soluzione esatta che è già costante a tratti. 
D'altra parte, se la soluzione del problema di Riemann contiene onde centrate di rarefazione o onde 
miste, allora approssimeremo un’onda di rarefazione (wf_j, wf) con un “ventaglio di rarefazione” 
costante a tratti, ottenuto inserendo degli stati intermedi tra wf_j e wf che corrispondono a piccole 
discontinuità che viaggeranno con una velocità vicina alla velocità caratteristica A;(wf_,). Natu- 
ralmente le discontinuità dei fronti d’onda del ventaglio di rarefazione non soddisfano le condizioni 
di (RH) esatte (in quanto gli stati intermedi non sono connessi da una curva di shock), ma soddis- 
fano delle condizioni di (RH) approssimate. La soluzione approssimata u così definita può essere 
prolungata fino al primo istante t in cui interagiscono (si incontrano) due (o più) fronti d’onda. A 
questo punto, poichè u(t1, -) è ancora una funzione costante a tratti, possiamo di nuovo risolvere 
in maniera approssimata i corrispondenti problemi di Riemann con una funzione costante a tratti, 
così come abbiamo fatto al tempo t = 0. La soluzione può essere quindi ulteriormente prolungata 


fino al successivo istante #2 in cui si verifica un’altra interazione tra fronti d’onda, etc ... 


Per sistemi generali di n leggi di conservazione il principale problema tecnico posto da questo 
procedimento è dovuto al fatto che il numero dei fronti d’onda può diventare infinito in tempo finito. 
Infatti ad ogni generico punto di interazione ci saranno in generale due fronti d’onda entranti, 
mentre il numero dei fronti d'onda uscenti sarà n (se tutte le onde elementari che risolvono il 
problema di Riemann sono shock o discontinuità di contatto) o maggiore di n (se la soluzione del 
problema di Riemann contiene qualche rarefazione od onda mista). Per superare questa difficoltà 
è necessario modificare l’algoritmo fin qui descritto adottando due modi differenti di risolvere il 
problema di Riemann: uno è quello appena richiamato che introduce vari fronti d’onda nuovi. 
L’altro invece è un algoritmo semplificato che determina una soluzione contenente un numero 
minimo di fronti d’onda. In questo secondo caso, tutte le onde nuove che verebbero create sono 
raggruppate insieme in un unico fronte d’onda artificiale (fronte d’onda non fisico) che viene 
fatto viaggiare con velocità maggiore di quella di tutti gli altri fronti d’onda. Questo algoritmo 
semplificato viene usato quando la grandezza dei fronti d’onda interagenti (la grandezza del salto) 
è più piccola di una certa soglia fissata. 


Utilizzando il procedimento fin qui descritto si possono costruire soluzioni approssimate di (1)- 
(2) definite per tutti i tempi t > 0. L'argomento principale della dimostrazione di questo risultato 
si basa sulla derivazione di stime a-priori della variazione totale ottenute introducendo un 
funzionale del second'ordine (rispetto alla grandezza dei fronti d’onda presenti nella soluzione) 
che misura l’interazione potenziale tra fronti d’onda. Argomenti di compattezza permettono 


poi di mostrare che una sottosuccessione converge ad una soluzione debole entropica di (1)-(2). 


Possiamo allora enunciare il risultato generale di esistenza per il problema di Cauchy (1)-(2) 


con condizione iniziale a variazione totale piccola. 
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Teorema 1. Sia il sistema (1) strettamente iperbolico e supponiamo che i campi caratteristici 
siano 0 (LD) o (GNL), oppure (NGNL). Allora esite una costante $ > 0 tale che, per ogni 
condizione iniziale © € L'(R; R"), con 


Var. Tot(u) < $, (46) 


il problema di Cauchy (1)-(2) ammette una soluzione debole entropica definita per tutti i 
tempit > 0. 


x 


Questo risultato è stato ottenuto per la prima volta, per sistemi con campi caratteristici 
(LD) o (GNL), usando lo schema di Glimm in [G]. Le idee principali sul tracciamento dei fronti 
d’onda sono state introdotte in [Da] per equazioni scalari, in [DP1] per sistemi di due leggi di 
conservazione, e poi estese al caso generale di n leggi di conservazione (sempre sotto le ipotesi che 
i campi caratteristici siano (LD) o (GNL)) in [BJ], [Bi], [Ri]. 


Il caso non convesso è stato estensivamente studiato da T.P. Liu che ha dapprima determinato 
la soluzione del problema di Riemann in [LI1], [LI1], e poi ha dimostrato l’esistenza globale di 
una soluzione del problema di Cauchy, usando lo schema di Glimm, in [LI4]. In [AM3] abbiamo 
introdotto un algoritmo di tracciamento dei fronti d’onda che permette di ottenere il risultato 
di esistenza del Teorema 1 per sistemi di due leggi di conservazione con campi caratteristici 
(NGNL). L'estensione di questo risultato al caso più generale di sistemi di n leggi di conservazione 
con campi caratteristici (NGNL) è trattata in [AM4]. 
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